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Abstract 



We determine the position and shape of the melting line in a layered super- 
conductor taking the electromagnetic coupling between layers into account. 
In the limit of vanishing Josephson coupling we obtain a new generic reentrant 
low-field melting line. Finite Josephson coupling pushes the melting line to 
higher temperatures and fields and a new line shape B m oc (1 — T/T c ) 3//2 is 
found. We construct the low-field phase diagram including melting and de- 
coupling lines and discuss various experiments in the light of our new results. 
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Since its proposal in 1988 JT| , vortex-lattice melting in bulk type II material has become 
a central topic in the phenomenology of high temperature superconductors. The order, po- 
sition, and shape of the transition have been investigated theoretically |2j as well as experi- 
mentally 0| by a large number of authors. Most recently, the main interest is concentrating 
on the phase diagram of the strongly layered B^S^CaiCh^Og (BiSCCO) superconductor 
which is being investigated by means of /iSR j4|, neutron scattering ||, SQUID magnetom- 
etry ||, and Hall-sensor arrays 0], probing the melting- and/or decoupling transition in 
these materials. It turns out that the most interesting regime is the low-field part of the 
phase diagram with B < 1 kG, where the electromagnetic interactions between the layers 
becomes relevant, and it is the purpose of this letter to derive and analyze the vortex-lattice 
melting transition in this regime, taking full account of electromagnetic coupling. 

The importance of electromagnetic interactions, contributing to the stiffness of individual 
vortex lines, has been realized before within the context of vortex-lattice melting in the dilute 
limit ||, where the transition line exhibits a reentrant behavior (lower branch of the melting 
line). As we will show below, the electromagnetic interaction also influences the behavior 
of the upper branch of the low-field melting line and even may change its shape from the 
usual B^(T) oc (1 — T/T c ) 2 behavior to a new power-law B^' 3 (T) oc (1 — T/T c ) 3 / 2 within a 
large part of the phase diagram — this is the new and central result of this paper. 

Our analysis below is based on the continuum elastic description of the vortex lattice 
combined with the Lindemann criterion, stating that the lattice will undergo a melting 
transition once the mean thermal displacement {u 2 ) l J^ becomes comparable to the lattice 
spacing a Q ~ (Qo/B) 1 / 2 , (u 2 )H 2 /a |T m ,s m ~ c L . The Lindemann number c L is usually chosen 
to be a constant of order c L ~ 0.1 — 0.3. Though not rigorous, the Lindemann-type melting 
scenario has proven very useful and reasonably accurate in predicting the positions of first- 
order melting transitions in general and the line shape of the vortex-lattice melting transition 
in particular. 

A well known limiting case, where strong fluctuations due to dimensional reduction drive 
a vortex-lattice melting transition, is the superconducting film (2D dislocation-mediated 
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Kosterlitz-Thouless melting scenario, see 0) and we will begin our analysis with this el- 
ementary building block of a layered superconductor. Next, we consider a layered system 
with electromagnetic coupling and derive the shape of the reentrant melting line in this 
limit. Finally, we account for the Josephson interaction between the layers producing a 
finite anisotropy parameter e 2 = m/M < 1, where m and M denote the effective in-plane 
and c-axis masses. Our results are illustrated in Fig. 1, where we show the shape of the 
vortex-lattice melting line as it evolves from the 2D isolated layer, to the electromagnetically 
coupled system of layers, to the Josephson coupled bulk anisotropic superconductor. 
Our main task is the calculation of the mean-squared thermal displacement |L0] 
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(in (|TJ) we neglect a second contribution to (u 2 ) t h involving lattice compression and keep 
only the main term). Here, e Q = ($ /47rA) 2 denotes the basic energy scale of the continuum 
elastic theory, $ Q = hc/2e is the flux quantum, A denotes the planar London penetration 
depth, and £ is the planar coherence length. The second term in the single vortex tilt 
c c AA is due to the electromagnetic coupling between the layers and is the only term in C44 
surviving the limit e — > (layer decoupling). The electromagnetic contribution to the tilt 



modulus is strongly dispersive and produces the large stiffness E\ ~ e /2 of the vortex lines 
in the long-wave-length limit k z < 1/A. With increasing k z the electromagnetic stiffness 
decays oc 1/X 2 k 2 and the line tension crosses over to the well known result ei ~ e 2 e for 
the anisotropic superconductor as k z increases beyond 1/eX (note that this residual tension 
is due to the Josephson coupling and is relevant only for eX > d, where d denotes the 
layer separation). The expression given in (f|) is valid for small displacements, in the elastic 
regime. For large displacements uk z > 1 the logarithm in the second term of (f|) should be 
cut on 2X/u rather than Xk z fT2j. In our analysis below we then replace the logarithm by 



the factor X 2 k 2 /(1 + (3X 2 k 2 ) with (3 = l/ln(l + AX 2 /c 2 a 2 ) producing a smooth interpolation 
between the hard and soft tilt modes at large and small wave-lengths, respectively. 

We start with the analysis of an individual layer (we use the definition X 2 jd = X 2 /d s with 
A s and d s the penetration depth and thickness of the superconducting layer). Dropping the 
tilt energy in ([]]), the integral over k z provides a factor 2it/d and cutting the .fT-integration 
on a few lattice spacings we obtain the ratio (u 2 ) t \ 1 /a 2 « T ^nc^da 2 . Within our simple 
Lindemann approach we then can reproduce the correct result « a^dc^j^K for the 
dislocation-mediated melting temperature if we choose a Lindemann number c L = 1/2^/2-k « 
0.1. The high-field part (a < A) of the melting line is field-independent, 

and using parameters typical for the layered high-T c superconductors, T c w 100 K, A 2 (T) w 
Ag/(1 - T 2 /T 2 ) with A w 1800 A, and d = 15 A, we obtain e q (T = 0)d w 10 3 K and 
« 20 K. The low-field part (a > A) of the melting line is dominated by the exponential 
decay of the shear modulus and we obtain the result 
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The result for the melting line of an isolated layer is illustrated in Fig. 1. Note that the 
result (|6|) is a somewhat artificial construct, as we have used the low-field expression in (g) 
for our analysis. In this way we illustrate the behavior of the melting line in the absence 



of any interlayer coupling (neither electromagnetic nor Josephson) while keeping the shear 
modulus of a translation invariant system along the field axis. On the other hand, the 
analysis of a 2D film involves a different shear modulus |13[] which softens only at very low 
fields (a Q > Acfr = 2\ 2 /d) following the power-law behavior cqq ~ 0A6e o X e s/a 3 oc B 3 ^ 2 
rather than the exponential behavior used in our analysis. 

Next we consider a finite electromagnetic coupling between the layers while keeping e = 
(no Josephson coupling). In the high- field regime (a < A) the shear term in ([!]) dominates 
over the tilt energy and we recover the field independent 2D-result (f|). For small fields with 
a Q > A the tilt energy becomes relevant and the Lindemann criterion reads 
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with 5 = 2cqqX 2 /Eo- Here, the first term originates from the soft tilt modes with k z X > 1, 
whereas the second term involves the long wave-length modes hardened by the electromag- 
netic coupling. This second term becomes relevant only at very small fields a / A>1, where 
the shear modulus is exponentially small, 5 oc exp(— a /\). The result (|?p provides a lower 
branch of the melting line which is limited by soft shear and hard tilt, 
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and no solid phase can exist at high temperatures beyond T x . Using typical parameters 
for the layered high-T c materials and adopting a value c L ~ 0.1 for the Lindemann number 
we find T x close to T c , 1 — T x /T c w 0.05 (in (^) we have introduced the 2D Ginzburg 
number G 2D ~ T c /e (T = 0)d ~ 0.1; the logarithms in (§) and (|H]) take typical values 



around 5-6). The reentrant melting line defined by (|j) and (||) is illustrated in Fig. 1: 
The electromagnetic coupling of the layers favors the solid phase and the low-field melting 
line develops the characteristic "nose-like" shape of a 3D system. Note that the point of 
reentrance ends up in the critical region close to T c . Since our approach accounts for the 
fluctuations of the phase-field via the thermal motion of vortices but neglects amplitude 
fluctuations of the order parameter our analysis breaks down in this regime. 

In the final step we account for the Josephson coupling between the layers producing 
a finite anisotropy parameter e > 0. This additional coupling becomes relevant whenever 
a Q ,X > d/e and again favors the solid phase, hence pushing the melting line further to- 
wards higher temperatures and fields. Evaluating the Lindemann criterion in the low-field 
regime (a > A) we recover the previous result (0) with the modification that the soft tilt 
modes are cut off on l/e\f]3\ instead of ir/d, leading to the replacement of ln(. . .)/47r£ by 
(dy/]3/7ieX)[\n(. . .)/Airf35 + 1] in the first term of (|7|). The lower branch of the melting line 
remains unaffected, whereas the upper branch of the low-field melting line takes the new 
form 

$ IXC 2 EE X / T 2 \ 3 / 2 



The crossing point of the lower and upper branch of the melting line is shifted towards 
higher temperatures, 
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For eXq < d the line 5J^ n,J goes over into the generic melting line B^ ,u as the temperature 
drops below T em w T c [l-p(EiiX /d) 2 } 1/2 . For the opposite case where eXo > d the generic line 
-B£f' u is completely hidden and i?^ n,J merges into the well known bulk anisotropic melting 
line as the field grows beyond $ /A 2 : At these fields the tilt energies are dominated 
by the dispersive bulk term c^ 4 ~ 4nE 2 E /a^K 2 (see Eq. (||D) and the Lindemann criterion 
provides the well known result 



Bi(T) « -^Aitcl-^- (13) 



At large fields (a < d/e < A ) the 2D result ([|) is recovered. 

The most interesting result is the new line shape B^ J ~ (1 - T/T c ) 3 / 2 , Eq. (fTJ), 
describing the low- field/high-temperature melting in a Josephson-coupled layered or highly 
anisotropic superconductor (small parameter e < d/X ). This new result is due to the 
electromagnetic coupling which dominates over the bulk-dispersive tilt modulus c^ 4 as well 
as over the single-vortex line tension e 2 e due to Josephson-coupling in this regime. The 
substitution of the old result (|13|) by the new expression (pT|) is particularly relevant in the 
strongly layered superconductors such as BiSCCO: The (1 — T/T c ) 3 / 2 power-law is valid 
provided that djn\f^e < A < a . Assuming e ~ 1/150, the second restriction implies 
T > 0.4 T c . In less anisotropic materials, such as YBCO with e ~ 1/5, this condition is 
much more stringent and the upper branch of the melting line is always described by the old 
result, Eq. (0). Note, however, that in YBCO the suppression of the order parameter close 
to the upper critical field H C2 becomes relevant and the melting line cannot be described in 
terms of a simple power law oc (1 — T/T c ) 2 any longer, see Ref. 0, Blatter and Ivlev, for a 
detailed discussion (in BiSCCO the melting line is far below H C2 and there is no suppression 
of the order parameter in this regime). 

It is instructive to compare the different low- field melting lines as given by Eqs. (|9|), (|TTD, 
and (|i~3l) . A quick inspection gives the ratio B™' n / B^ = [dj Xec L ) 2 T/S7ij3e d, which is of 
order unity taking the above parameters for BiSCCO and using e = 1/150, a value often 
quoted in the literature |,§. Similarly, B^ 3 /B 3 m = (d/ Xec L )T / '16 y/pc h e d w a[T 2 /T c {T c - 
T)] 1 / 2 , were again a ~ 1 if we use the above parameters for BiSCCO. The comparison of 
experimental data for the irreversibility or melting line with the theoretical prediction is 
often used to extract an estimate of the anisotropy parameter e, particularly in the strongly 
layered materials Following up the above discussion we draw attention to an important 
problem with this procedure: If the anisotropy parameter is very small, say e < 1/500, the 
(upper branch of the) low-field melting line (where the comparison theory/experiment is 
carried out) is dominated by the electromagnetic coupling and no anisotropy parameter 
can be extracted. The analysis of the melting line can provide a reliable estimate for the 
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anisotropy parameter only if e is large enough such that either the bulk result is valid 
or the mixed electromagnetic/ Josephson result ( |TTD can be identified via its particular line 
shape. 

It is generally believed that the low-field melting transition takes the vortex lattice into 
a liquid of vortex lines. It then has been proposed that this line-liquid transforms into a 
pancake-liquid in a second transition where the layers decouple, see Refs. |TT1JT0|] . A simple 



estimate for the position and shape of this decoupling line is obtained in the following way: 
thermal wandering of the vortex line over a distance L produces a displacement amplitude 
(u 2 )th — LT/ei. The layers decouple when the mean thermal displacement between line 
segments in neighboring layers becomes of the order of the lattice spacing, (u 2 )l(^(L = 
d) ~ a . For a Josephson-coupled system the line tension is Si ~ e 2 e Q and we obtain the 
well known result B 3 Ac ~ \<& /(d/e) 2 ] e d/T oc 1 — T/T c . However, for small anisotropy the 
electromagnetic coupling dominates at low fields and using the short wave-length elasticity 
Si ~ e (d/X) 2 we find that the decoupling line follows the melting line B^' n . We then 
obtain a phase diagram where the decoupling and melting lines are separate transitions at 
low (T < T™) and high (T > T cm ) temperatures but close up in between. 

Recently, a first-order phase transition has been observed in the low-field regime of a 
strongly layered BiSCCO superconductor f7j. The jump in the magnetization observed at the 
transition can be associated either with a vortex-lattice melting- or with a layer-decoupling 
transition. Fits using a (1 — T/T c ) 155 (melting) or a (T c /T — 1) (decoupling) power-law 
behavior produce a satisfactory agreement with the data over most of the measured tem- 
perature interval f?|. Our new result ( |TTD then is in good agreement with the measured 
power-law behavior based on the melting scenario. Whether the observed transition indeed 
can be attributed to a first-order melting transition remains to be shown, however. 

In conclusion, we have determined the position and shape of the melting line in a lay- 
ered superconductor taking the electromagnetic interaction between the layers into account. 
Whereas the electromagnetic coupling is irrelevant at fields B > $ D /A 2 , new results for 
the melting line have been obtained in the low-field regime B < $ Q /A 2 . In this regime, 
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the electromagnetic coupling produces a stiffening of the vortex line at long wave lengths 
k z < 1/A. Both, the lower and the upper branch of the reentrant melting line are effected 
by this stiffening and a characteristic "nose"-shaped 3D melting line is found even in the 
absence of Josephson coupling between the layers. Accounting for an additional Josephson 
coupling, the upper branch of the melting line is pushed out to higher temperatures and 
fields and takes on a new characteristic line shape oc (1 — T/T c ) 3//2 , as observed recently in 
a BiSCCO superconductor @. The results are crucial for an accurate understanding of the 
low-field phase diagram of layered superconductors. 

We thank E. H. Brandt for helpful discussions and the Swiss National Foundation for 
financial support. 
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FIGURES 

FIG. 1. Low- field phase diagram of a strongly layered superconductor. Reduced units 
b = 5/($ /Aq) and t = T/T c have been used and the Lindemann number c L = 0.1 has been 
chosen (else parameters appropriate for BiSCCO have been adopted, see text). The thick lines 
show the results for the isolated 2D layer and for the electromagnetically coupled system. The 
thin lines incorporate the effect of a finite Josephson coupling between the layers for anisotropy 
parameters e = 1/500, 1/150, and 1/50. The inset shows the same results on a logarithmic field 
scale, where the reentrant behavior of the melting line becomes more visible. The dotted line traces 
b(t) = B/l<t> /\Ht)]. 
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